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Abstract 

We study the reduced dynamics of a spin (qubit) coupled to a spin-boson environ- 
ment in the case of pure dephasing. We derive formal exact expressions which can 
be cast in terms of exact integro-differential master equations. We present results 
for a SB environment with ohmic dissipation at finite temperatures. For the special 
value of the ohmic damping strength K = 1/2 the reduced dynamics is found in 
analytic form. For K <C 1 we discuss the possibility of modulating the effect of 
the SB environment on the qubit. In particular we study the effect of the crossover 
to a slow environment dynamics, which may be triggered by changing both the 
temperature and the system-environment coupling. 
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1 Introduction 



The possibility of controlling the dynamics of a complex quantum system 
would open a wide scenario for both fundamental and applied physics. New 
physical properties emerging at the nanoscale may yield new paradigms for 
nanoelectronics, an example being quantum information processing [1]. The 
basic operations of a quantum logic device can be described in its compu- 
tational space defined by a set of degrees of freedom which can be prepared 
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(write) and measured (read) and of course evolve coherently. However the 
Hilbert space of the device is much larger than the computational space, and 
effects of the additional degrees of freedom may be relevant for the reduced 
dynamics we aim to control [2,3]. These environmental degrees of freedom 
may be "internal" to the device, which is a many-body object, or "external", 
belonging to the read/write circuitry or to auxiliary devices (buses or other 
idle or working qubits). They determine loss of coherence, which is ultimately 
due to the entanglement between system and environment [2]. 



Much of the effort in the analysis of decoherence in quantum logic devices has 
been focused on model environments made of many weakly coupled harmonic 
oscillators [2,4]. The key feature of these linear models is that the effects on the 
reduced dynamics of the system can be described in terms of simple statistical 
information on the environment, spectral densities or power spectra, even if a 
detailed knowledge of its microscopic structure is not known [5,6,7]. However 
a slightly more detailed knowledge of the environment dynamics is needed in 
certain situations, for instance if the environment has memory on the time 
scale of the system dynamics [8,9]. 



Bistable impurities located close to a solid state device which produce 1/f 
noise [10] are a physical realization of an environment with memory. They have 
been recognized as the major source of decoherence in recent experiments [11] 
demonstrating quantum coherent dynamics in superconducting devices. Such 
an environment has been modeled in Ref. [9], with a set of Fano impurities [12] 
which randomly polarize the device. A peculiar effect of this environment, 
related to its non linear nature, is that the reduced qubit dynamics shows 
various qualitatively different regimes depending on the switching rate of the 
impurities. In particular dephasing due to slow impurities is not described by 
a single rate and memory effects of the environment enter the qubit dynamics. 



The main goal of this work is to discuss a situation where the reduced dy- 
namics of the system is able to probe different regimes of the environment 
dynamics. To this end we introduce a model which exhibits tunable crossover 
effects. We take as the environment an impurity described by a spin-boson 
(SB) model [5,6]. The Hamiltonian is introduced in section 2, where we also 
discuss applications to physical problems. In section 3 the formal solution for 
the off-diagonal element poi(t) of the RDM is found, in the limit of pure de- 
phasing and ohmic damping in the scaling regime. In section 4 a generalized 
master equation for poi (t) is derived and the analytic structure of the kernels is 
exploited in section 5 to study the character of the behavior of poi(0- Finally 
in section 6 we sketch a phase diagram for this problem. 
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2 Model 



We study the following Hamiltonian (h = 1) 



H = --(J Z —cr x - - r z a z + H SB , (L) 

H S b = ^ r z - ^ t x - ]-t z X + w a oi°a 5 X = J2^a(ai + a a ) . 

1 1 1 a a 

It describes two spins (<x and f) coupled to each other by the interaction v. 
The spin f is also coupled to a collective coordinate X of a set of otherwise 
noninteracting harmonic oscillators (operators a) a and a a ). The Hamiltonian 
Hsb represents a spin-boson model [5,6], which describes a two state system 
with tunneling amplitude A and bias e$B coupled to a set of harmonic oscilla- 
tors. This set of bosons is characterized by the spectral density of the operator 
X, given by G(u) = J2 a ^^(^ — w„). The Hamiltonian Eq.(l) is very rich and 
may model many different physical problems. In this work we will regard a 
as describing a qubit coupled with an environment consisting of a bistable 
impurity (r) described by the SB model, Hsb- We are primarily interested in 
the reduced dynamics of the qubit, resulting from tracing out all the environ- 
mental degrees of freedom, i.e. the harmonic modes and the spin r. Physically 
a may represent a Josephson quantum logic device where external voltage 
gates and magnetic fields are used to tune the effective charging energy e and 
Josephson coupling Ej [4,13]. The environment may be a charge close to the 
device undergoing phonon-assisted tunneling between two positions. We will 
consider an ohmic spectral density G(u) [5], i.e. we choose a set of free bosons 
such that the power spectrum of the collective coordinate X is given by 



oo 

Sx(w)= I dt-(X(t)X(0)+X(0)X(t)) 

— oo 



= 7T G(u) coth ^ = 2tt K M e^l/" c coth ^ , (2) 

where u c represents the high frequency cut-off of the harmonic modes. The 
case of ohmic spectral density is interesting because the SB model has a rich 
phase diagram resulting in an environment where memory effects may be 
tuned by the temperature T. We will present results for K <^ 1 focusing on 
nonzero temperatures, T > A r — A(A/cu c ) K ^ 1 ~ K \ where A r is the renormal- 
ized tunneling amplitude in Hsb [5]. Physically this regime describes a SB 
impurity active in producing decoherence for the qubit a. From a technical 
point of view this regime requires some care. For instance the dephasing rate 
for (t z ) using the standard secular master equation [8,5] in the couplings \ a 
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(for e S B =v = 0) reads Sx(A r )/4 [5]. For T > A r it reduces to 
7 (T) = nKk B T, 



(3) 



this result being valid only if 7(T) <C A r [8,5]. Thus for large enough temper- 
atures an analysis which goes beyond the standard secular master equation is 
needed [16,5]. 

We mention finally other physical problems where the model we study may 
be relevant. First of all the Hamiltonian Eq.(l) may describe a system of 
two-qubits (a and f) subject to gaussian noise, the ohmic case describing 
the effect of a resistive circuitry. The low temperature regime, which may 
describe a working two-qubit gate, has been investigated using the standard 
secular master equation [14]. Instead for temperatures T > A r our model 
may describe a working qubit (a) coupled with a idle qubit (f) not perfectly 
switched off. Finally the spin a may represent a measuring device [15] for the 
mesoscopic system described by Hsb- This picture is an alternative framework 
for our work, namely the reduced dynamics of a provides the result of the 
measurement on a SB system. If the coupling v is weak and the dynamics of 
the SB is fast enough only the power spectrum of the coupling operator vr z 

oo 2 

StM = Jdt V - ((r z (t)r z (0) + r z (0)r z (t)) - (r,) 2 J (4) 

— oo 

enters the dynamics of a, whereas if the SB has a slow dynamics the spin a is 
able to detect also details of the dynamics of r which go beyond S T (u). 



3 Formal solution for the reduced density matrix 

We consider the pure dephasing regime Ej = where [H, a z \ = 0, so no relax- 
ation process takes place and dephasing is due to fluctuations of the energy 
difference of system, i.e. a Z) eigenstates. Hence the diagonal elements of the 
reduced density matrix of the qubit in the a z basis (populations) are constant, 
whereas off-diagonal elements (coherences) decay. All the information on the 
reduced dynamics is included in the coherences. In this section we derive an 
exact formal expression for the coherence pm(t) = ( | p(t) | 1 ) . 

We choose a factorized initial density matrix W(0) = p(0) <g> wsb(Q), with 
the spin-boson environment in the general state wsb(0). Product initial states 
are relevant in quantum computing, one of the main requirements being the 
possibility to prepare the qubit in well defined initial states. The coherence 
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Poi(t) is thus related to the following correlation function involving the spin- 
boson variables 

poi(f) =Poi(0)e-'Tr SB {e- lH -^ SB (0)e^} = p O i(0) e iet C-+(t) (5) 

where H± = H S b ± |t z . We evaluate the correlator C- + (t) using a functional 
integral approach. In the derivation we consider a dimensionless continuous 
"coordinate" q, and specify later to the spin states (r). The trace over the 
environmental degrees of freedom in the coordinate basis | q) \ {x a } ) reads 

C_ + {t)= j dq f J dq4q\ Jdxfdx i dx' i K_(q f ,x f ,t;qi,Xi,0) 

( qi, Xi | w SB (0) | ql, x\ ) K\{q s , x f , t; q\, xj, 0) (6) 
where we have introduced the compact notation 

x = {x a } , dx = Y[ dx a . (7) 

J a J 

The Feynman propagators K± are associated to the two Hamiltonians H± and 
can be expressed as a real-time path integral 

i(S ± [q] + S B [x] + S SB [q,x})] , 

with the constraints q(0) = qi, q(t) = qf and x(0) = Xi, x(t) = Xf . S B [x] 
and S SB [q, x] are the actions corresponding to the free boson and the spin- 
boson interaction terms in (1) respectively. We now specify to an initial non 
equilibrium state of the factorized form u>s_b(0) = w T (0) <g> wp, where the 
harmonic oscillators start from the thermal equilibrium state wp, and w T (0) 
is a general density matrix for the system of coordinate q. Then C-+(t) reads 



K±(qf,Xf,t;qi,Xi,0) = DqDxexp 



C-+{t) = J dq f J dqidq'i ( qi | w T (0) \ q[ ) J(q f , q f , t; qi, q\, 0) (8) 
J(qf,q f ,t;q t ,q' l ,0)= j ' DqDq' exp [i(S_[q] - S + [q'])} F FV [q,q'\ (9) 

with the path integration constraints q(0) = qi, q(t) = qf, q'(0) = q[, q'(t) = 
qf . The influence of the harmonic modes enters the generalized "propagating 
function" J(q/,qf, t; qi, q\, 0) via the Feynman- Vernon influence functional 



t t' 



FFvM=vv{j dt'J dt" W)Q\t' - t")i(t") 

+ ii{t')Q"{t'-t")r ] {t")\}, (10) 



o o 



5 



where we have introduced the linear combinations 



V(t) = [q(t)+q'(t)}/2, m = l<l(t)-q'(t)]/2, 

describing respectively propagation along the diagonal of the density matrix 
and off-diagonal excursions. For a spin these paths are piecewise constant, 
the time intervals in which the spin is in a diagonal state 77 = ±1 are called 
sojourns, while the time intervals spent in an off-diagonal state £ = ±1 are 
referred to as blips. Depending on the initial condition for the spin, the initial 
state can be either a sojourn or a blip. We now choose for the spin a diagonal 
initial state in the "coordinate", i.e. r 2 , basis w T (0) — \ l + ^5p(0)r z , where 
Sp(0) denotes the initial population difference between the r z eigenstates. Then 

all paths that contribute to C h (t) in Eq.(8) start out and end up in a sojourn 

state, q(0) = q'(0) = q { -> 77(0) = 77 = q h q(t) = q'(t) = q f -> 77(f) = rf = q f . 
Thus a general path is characterized by 2m transitions at the intermediate 
times tj (tj — 1, 2, . . . , 2m) and it is parametrized by 

m 

i {2m) (t>) = £ & [e(f - ^-0 - e(f - t 2j )\ , 

j Z' (11) 

r 1 {2m Kt') = Y.mm'-^)-e{t'-t 2]+l )] , 

j=0 



where t = and t2m+i = £■ Inserting the path (11) into S± in Eq.(9) we get 



exp 



1(5- [g] - 5+ [?'])] = Wi' (-—) B m expliv^VjSj 



12 



where B m = exp —iesB Y^Li ij r j is the usual bias factor which depends only 
on the blip lengths Tj. Notice that due to the two propagators K± which occur 
in Eq.(6), we get dynamic contributions which depend on the sojourn lengths 
Sj in the spin actions (last term in Eq.(12)). 

Substituting the paths (11) into Eq.(10), the influence functional reduces to 



(13) 



The factor G m contains all the inter- and intra-blip-interactions, 



G m - exp [ - J2]Ll Q'2j,2j-1 - T,™=2 Efc=l £j6Aj,fc] , 



A-j,k — Q'lj,2k-\ + Q'2 



!2j-l,2k 



Q'2j,2k Q'l 



2j"-l,2fc-l ) 
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the sojourn-blip interactions are captured by the phase factor H m given by 

m j—1 

H m = exp ZjVk X J,k , ] 

j=l k=0 

Xj,k = Q2j,2k+1 + Q2j-l,2k ~ Q2j,2k ~ Q2j-l,2k+l ■ 



The kernels Q' and Q" are the real and imaginary parts of the second integral 
of the bath correlation function 



oo 



Q(t) = Jdu {coth(^) (l - cos(cjt)) + isin(cut)} . 

o 

and we have used the notation Q h ^ = Q{t h — tj). 

We proceed considering the important case of an ohmic bath Eq. (2), where 
the correlator Q(t) in the scaling limit u c t 3> 1, reads 

Q(t) = 2K In sinh ^ + i 7rKsgn(t) . (14) 

Using the above scaling form each sojourn interacts only with the subsequent 
blip, so that H m simplifies to 

TO— 1 

H m = exp [iirK J2 Vktk+i] • (15) 

A:=0 

Taking into account this features, the sum over histories of paths of Eq. (9) is 
represented by 

(i) the sum over all possible numbers of steps, 

(ii) the time-ordered integrations over the corresponding flip times {tj}, 
(in) the sum over all possible arrangements of blips {£j = ±1} and sojourns 

{%-+!}• 
The correlator C_ + (i) becomes 

C-+ (t) = E w v e** + E V E rf x 



E 



M -— J / V 2m {tj} Y G m B m Y H m exp [iv ^ , (16) 

where = ( 77 1 u> T (0) | rj ) and denotes the sum over all paths rjj subject 

to the constraints r] = i], i] m = r]'. The first term in Eq.(16) comes from the 
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path without any transition between the initial and final sojourns, and occurs 
for i] = i]'. The sum over {rjj}' can be easily performed using the scaling form 
(14) for H m , this leads to 



Y, H m exp [iv Y rjjSj] = e mKr ^ e lv ^ So+v ' Sm) 2 m ~ l D m (17) 

{%}' 3=0 

where D m = YYJL2 cos ( n K£j +vsj-i). Whereas for v = the sum over the 
intermediate sojourns provides a contribution which only depends on the blip 
states £j, because of the coupling the sum over rjj gives also a dynamic term 
dependent on the sojourn lengths Sj. Thus we finally get 



C-+(t) = E w v eiVVt + A (t) + iB(t) (18) 
v 

A(t) = Wvw v Yv'A(v,v',t) (19) 
z v v' 

B(t) = W^J2v'B( V , V ',t) (20) 
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Z T) T)' 



where 



00 / A 2 \ m \ 

A(ri, 1/, f) = £ -— / V 2m {t 3 }e^ s ^'^ £ G m B m D m cos(ttK) (21) 

m=l V ^ / jj {<-} 
oo / A 2 \ m ' 



m=l \ / o {£,} 

Eq.(5) with Eqs.(18) - (22) is the exact formal expression for the coherence 
Poi(t), valid for Ohmic damping in the scaling regime. 



4 Generalized Master Equation 



For quantum dissipative systems it is generally useful to formulate the dy- 
namics of the open system by means of appropriate master equations for the 
reduced density matrix. For the spin-boson problem generalized master equa- 
tions relating the averages of the spin components r a have been obtained in 
Refs. [5,16]. In this section we show that the series expressions for A(r),r)',t) 
and B(r),r)' ,t) given in Eqs. (21,22) can also be cast into the form of exact 
integro-differential equations. 
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To this end we have to find the kernels, i.e. the irreducible components in 
the formal expressions (21,22). Irreducibility of a kernel means that it cannot 
be cut into uncorrelated pieces at an intermediate sojourn without removing 
correlations due to the harmonic bath, across this sojourn. Following the lines 
of Refs.[5,16] we define irreducible influence functions subtracting from the 
inter- and intra-blip interaction factors G m all the reducible components 



G n — G n ^ ' ( l)"' ^ ' G mi . . . G m .5 mi j rm2 j r _^ m .^ n , 

j=2 m\...mj 

where the inner sum is over all positive integers rrij. By inspection of (21,22) 
we find that the irreducible kernels are 



K l {t-t')=K, l {t-t') 

(TrK)f^(~\ j 'dtto-i- j ' dt 2 J2G n B n D n (23) 



2 \ n t 



cos( 

n=2 \ / ,, ,, 



K 2 {t-t')=K 2 {t-t') 

+ sin(TrK) p ( ~ ) j dt 2n _ x . . . / df 2 £ iiG n B n D n , (24) 



V 



where the lowest order contributions do not depend on the coupling v and 
coincide with the lowest order kernels entering the spin dynamics in the un- 
coupled case [5,16] 

/Ci(f - = A 2 cos(7r^)G'i(t - t') cos\e SB (t - t')\ (25) 
K, 2 {t - t') = A 2 sin(7rX)G 1 (i - t') sm[e SB (t - t')] . (26) 

With the above the generalized master equations for A(r/, r/', t) and B(r/, rf , t) 
read 



d 



—A = irj'vA - J dt' [e vr ' vt ' + A] K^t - t r ) 
o 

t f 

+ ivj dt'e-^W J dt"[ V 'K l {t' - t") - K 2 {t' - t") }A(t") (27) 

o o 

t 

B = ir]'vB - J dt'lK^t - t')B - iK 2 {t - t')e^ vt '} 
o 

t f 

+ ivj dt'e-^W J dt"[r ] 'K 1 {t' - t") - K 2 (t' - t") }B(t") . (28) 



d_ 

dt 
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In the absence of coupling, v — 0, A(f],i]',t) and B(r],r]',t) are simply related 
to the components of the average (r z ) t symmetric and antisymmetric under 
bias inversion e SB -> ~^sb, {r z )t = P s (t) + P a (t), via A (r),r)',t) = P s (t) - 1 
and B (rf,rf\t) = —iP a (t). In this limit Eqs. (27,28), reduce to the master 
equations satisfied by P s (t) and P a (t) [5,16]. 

Being of convolutive forms the integro-differential equations Eqs. (27,28) are 
conveniently solved in Laplace space. Performing the sums in Eqs. (19,20) the 
Laplace transform of the correlator C-+(t) reads 

C_ + (A) = -^ [X + K^-ivSpiO)] (29) 
D(X) = X 2 + v 2 + XK 1 {X)+ivK 2 {X), (30) 

where Ki(X), K 2 (X) are the Laplace transforms of the kernels (23,24) and 
Sp(0) comes from the initial condition chosen for the spin f. 



5 Reduced dynamics at finite temperatures 



Solving the pole equation (30) represents a formidable task, so in the follow- 
ing we shall focus on special damping and temperature regimes which are 
particularly relevant for quantum computing. We shall analyze dephasing at 
finite temperatures at two special coupling conditions: weak damping K <C 1 
and K — 1/2. In the second case the spin-boson model can be mapped onto a 
Fano Anderson model [12] which is suitable to study dephasing due to bistable 
fluctuators producing 1/ / noise [9] . 

We shall adopt an approximate treatment widely used in the spin-boson liter- 
ature [5,6], the noninteracting-blip approximation (NIB A). In the NIBA, the 
blip-blip interactions Aj^ are neglected, and also the sojourn-blip interactions 
Xj t k are disregarded except those of neighbors, k — j — 1, and they are approx- 
imated by Xjj-i = Q"{Tj). Thus the influence functions (13) factorize into 
the individual blip influence factors, and the irreducible influence functions 
(23) vanish for n > 1. For Ohmic damping in the scaling regime, each blip 
interacts with the subsequent sojourn only (see Eq.(15)), thus the only ap- 
proximation involved in the NIBA is the neglect of the blip-blip interactions. 
This approximation is justified in the regimes of our interest, see ref. [5]. 

Within the NIBA the series expression for the irreducible kernels (23,24) re- 
duce to the lowest order contributions (25,26). Evaluation of Eqs. (27,28), with 
these kernels give the reduced dynamics of the qubit in the NIBA. To make 
things simpler we shall consider the unbiased case esB = 0. 
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5.1 The case K =1/2 



For K = 1/2 the Laplace transform of the kernel Ki(t — t') takes the form 

oo 2 

k^i/2J 1K ' 2 u c ' ' 

o 

the zero coming from cos(7ri^) for .K" — > 1/2 being compensated from the short 
distance singularity of the intra-blip interaction lim^o e~ Q '^ « (uo c t) 2K . The 
pole equation (30) is readily solved and C_+(i) reads 



C_+(i) = Ae-^ (1 - a)7F * + (l-A)e _ s (1+a)7li,t . 



(31) 



This form is suitable for illustrative purposes since a and A have a very simple 
expression [17] 



a 



/2v\2 



1 r 2w 
A = — 1 + a - i — 5p(0) 
2a L 7f ■ 



(32) 



The behavior of C_+(i) has been analyzed in Ref. [9] where two qualitatively 
different regimes for the problem were identified. 

Fast environment- If 2v/jf <C 1 only one of the exponentials appearing in 
Eq. (31) dominates the long-time behavior and C-+(t) shows a single rate, 
given by |(1 — a) 7jt ~ v 2 / r yF, the standard Golden Rule result. The decay of 
C_+(i) at all times, described by T(t) = — ln{C_ + (t)}, is well approximated by 
simulating the SB model with suitable set of harmonic oscillators reproducing 
the power spectrum S T (u). In this case the decay is described by T(t) ps T osc (t) 
where [3,5] 

oo 

t-, /x f du _ . s 1 — cos cut 

r osc (t) = / — S T {u) - 2 . 33 





Here the power spectrum Eq.(4) is a Lorentzian S t {uj) = 2v 2r yp/ {j f + oj 2 ), due 
to the pure relaxation behavior of (t z ) for v — [5]. Thus in this regime the 
spin a sees the SB environment as an effective set of harmonic oscillators, and 
no detailed information or memory effects of the environment enters C ^(t). 

Slow environment- If 2v/ r yp ^> 1, C h (t) is determined by both the expo- 
nential terms in Eq. (31). Deviations from the weak coupling result T osc (t) 
are observed, in particular the short times behavior t < 7^' shows memory 
effects whereas for long times the coupling v does not enter the decay rates, 
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given by = 7f/2, but determines the frequency shifts, which reduce to ±v 
(see Fig. 5.1. a). 




Fig. 1. (a) The real part (solid lines) and the imaginary part (dots) of the rates 
entering Eq. (31) are plotted for finite temperature ksT ^> esB > 0. (b) The full 
expression for T(t) for various 2v/^f- deviations from the weak coupling result 
r osc (t) (dashed line) are apparent for < 2v. Here 5p = tgh [€sb/(2/cbT)] = 0.08. 



We mention that for K = 1/2 the same structure Eq. (31) is also found in 
more general situations, €sb,Ej 7^ [9,17]. Formulas for A and a are more 
complicated, for instance TZe a and 1m a are both non vanishing in the entire 
2v/~fF range and a third regime 2v/jf ~ 1 appears (see Fig. 5.1). However the 
analysis using the expressions (32) remains still valid, and in the following we 
will identify the different regimes for if < 1 from a similar analysis of the 
poles of C*_ + (A). 



5.2 Weak ohmic damping K <C 1 and finite temperatures 

For weak ohmic damping, K ^ 1 and temperatures hsT^ < ksT u c , where 
ksTb = \J + v 2 , the characteristic fluctuations of the harmonic oscillators 
in Hsb have no memory and Q(t) assumes the Markov form 

Q(t) = 2K{^- + In (^) } + Mr* sgn(f) . 

Using this latter expression the interblip interactions cancel out exactly, thus 
justifying the NIBA. 

In the absence of bias, = only the kernel /Ci(A) is non vanishing, and 
in the Markov limit it reads /Ci(A) = A^/[A + 27(T)], where 7(T) is given by 
Eq. (3) and we have introduced the temperature dependent tunneling term 
A T = A r (2nkBT / A r ) K . The correlator C-+(\) is readily found as 
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r rM _ [A + 2 7 (T)] [A_-njg(0)] + A| 

C - +(A) " A(A) (34) 

N(X) = (X 2 + v 2 ) [X + 2 1 (T)] + A 2 T X 

The qualitative behavior of C_+(i) is determined by the three poles found by 
solving the equation N(X) = 0. 

Low temperatures - At temperatures T5 < T < T the poles, in lowest order 
in T/T , are given by 



A » = "^ ! A V, = -^±V« 2 + ^. (35) 

The limiting temperature is given by /cbT p=j A t /(ttK) — v 2 /(2nKA T ) for 
w/A T < 1 and by k B T = ^ A 2 T + v 2 /(2nK) for w/A T > 1. The behavior of 
C- + (t) depends essentially on v/ A T . If v / A T <C 1 the reduced dynamics shows 
a single rate, corresponding to the dominant pole Ao ~ — 2(v/A<r) 2 7 (T). In- 
stead in the opposite regime, v/ At ^> 1, the two poles A1/2 ~ — (At/v ) 2 7(T)± 
iv dominate. The frequency shifts are given by ±t> as in the slow environment 
case of section 5.1. Thus for v/A T ~ 1 a crossover is expected in the behavior 
of C_+(t). 

High temperatures - If T ^> T , the leading behavior is given by 



A = -2 7 (T) ; A 1/2 = --4|r± 1 



A 4 

A ^ _ 4u 2 . (36) 



4 7 (T) 2\j4 7 2(T) 

In this case there is an intermediate temperature regime, A^Tq <C /c^T < 
ksT + = A T /(4nKv), where we may have three real roots. In this regime, 
which occurs only if v/A T < 1, the dominant contribution to C_ + (i) comes 
again from a single pole, Ai ps — 2(v / A T ) 2 l{T). Instead for temperatures 
T > T + the two complex conjugate A1/2 are the dominant poles. The decay 
rate displays the Kondo behavior A^/[4 7 (T)] oc T 2K ~ l and does not depend 
on v. The coupling t> enters the frequency shifts, which approach ±t> if the 
temperature is further increased. Thus in the region T»T and v < A T we 
find a crossover which may be triggered by both v and T, the crossover line 
being given by 

A 2 

v . (37) 



Finally for larger values of v > At the system stays in the regime where 
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the dominant rates are the two complex conjugate A1/2 ~ —Aj./[A'y(T)] ± iv. 
Again v does not enter the decay rate but determines the shifts ~ ±v. 



6 Discussion of the results and conclusions 

The analysis of the poles of C |-(A), Eq.(34) allows to extend the qualitative 

conclusion for K = 1/2 to the regime K <C 1. In this latter case we find that 
the dynamics of a shows a behavior which reflects the crossover of the SB 
environment to a slow dynamics regime. For K <C 1 the SB environment can 
be modulated by acting on both the coupling v and the temperature. Instead 
for K = 1/2 and esB — the crossover is triggered by the parameter v/'-fp only. 
This is simply due to the fact that for the relatively strong coupling K = 1/2, 
t z has a relaxation dynamics at any temperature. Thus at any temperature 
we may identify the correlation time of the SB environment r c pa 7^. The 
crossover to the slow environment regime occurs when 2v r c ~ 1 which marks 
the limit where the standard master equation approach in the coupling v would 
break down [8]. 

To understand the behavior for K <C 1 we focus on t) « A T and we iden- 
tify the correlation time r c of the SB environment from the dynamics of 
the SB model for v — 0. In this case r 2 displays a crossover from coher- 
ent to incoherent dynamics at a temperature kBT*(K) pa A r /nK. Below 
T* coherent oscillations at frequency ~ A T are slowly damped with rate 
given by 7(T) <C A T and we may identify r c ~ A T X . Thus the SB envi- 
ronment is fast, since 2vr c pa 2v/A T <C 1 and pmit) shows the single rate 
A ~ —2(v/Arp) 2/ ~f(T) pa S T (0)/2. If the temperature increases fairly above 
T* the dynamics of f is incoherent and relaxation occurs, with a slow rate 
7r(T) pa A^/27 (T) [5] which is identified with r c _1 . Slow environment effects 
are then expected to appear in the reduced dynamics of a at 2v/^ r {T) ~ 1 
which is indeed the crossover line Eq.(37). 

Finally we point out possible extensions of our analysis. A SB environment 
with subohmic (superohmic) spectral densities [5] could be treated along the 
same lines, but since the uncoupled SB dynamics is always slow (fast) we 
do not expect tunable crossover lines to appear. Our analysis can be easily 
extended to an environment made of an arbitrary number of SB impurities. 
In this case C- + (t) = where C^ + (t) is the individual contribution 

of the i— th SB impurity. This allows to model a simple environment which 
produces 1/ / noise by choosing a suitable distribution [10] of the parameters 
of the set SB impurities. For instance in Ref. [9], which focuses on K — 1/2, 
the standard distribution oc I/7F has been chosen for the relaxation rate of 
the SB. We notice that for K = 1/2 temperature has no effect in the regime we 
consider. If we allow for impurities with Esb 7^ the temperature determines at 
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most how many impurities in the set are active. Instead for a 1/ f environment 
modeled by a set of ohmic SB impurities with K <C 1 our results show a non 
trivial behavior as a function of the temperature, which may be relevant for 
recent observations [18] of temperature-dependent dephasing due to 1/ f noise. 
Finally if Ej ^ 0, a z is not conserved and relaxation occurs. The relaxation 
rate increases monotonically with Ej/e. We believe that the physical picture 
emerging from this work still applies to dephasing (see Ref. [17] for K = 1/2) 
at least for times smaller than the relaxation time. 
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